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4|k Based on Observation with the ESA Hipparcos satel- 
lite. 

Based on observations with the NASA/ESA Hubble 
Space Telescope. 

Abstract. Accuracy in the absolute position in the sky 
is one of the limiting factors for pulsar timing, and tim- 
ing parameters have a direct impact on the understanding 
of the physics of Isolated Neutron Stars (INS). We re- 
port here on a high-accuracy measurement of the optical 
position of Geminga (ruy = 25.5), the only known radio- 
quiet INS. The procedure combines the Hipparcos and Ty- 
cho catalogues, ground-based astrometric data, and Hub- 
ble Space Telescope (HST) Wide Field Planetary Camera 
(WFPC2) images, to yield Geminga's absolute position to 
within ~ 40 mas (per coordinate). Such a positional ac- 
curacy, unprecedented for the optical position of a pulsar 
or an object this faint, is needed to combine in phase 7- 
ray photons collected over more than 20 years, i.e. over 
2.5 billions of star' revolutions. Although quite a difficult 
task, this is the only way to improve our knowledge of the 
timing parameters of this radio silent INS. 
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1. Introduction. 

Geminga is the only known example of a bona fide pul- 
sar that cannot rely on well established radio timing tech- 
niques. The radio silence of the source implies the impossi- 
bility to measure its position and distance as inferred from 
the optimization of the radio timing parameters and from 
the dispersion measure of the radio pulses. To measure 
position, distance and timing parameters of Geminga, all 



other branches of astronomy had to be exploited in a 20 y 
long chase, recently summarized in Bignami and Caraveo 
(1996 and references therein). Briefly stated, the source 
was discovered in high energy 7-rays by the SAS-2 satel- 
lite in 1972, and studied in more detail by the COS-B 
mission. An X-ray counterpart has been proposed in 1983 
and an optical one, refining the position, in 1987/88. How- 
ever the breaktrough came with the discovery of the 237 
msec periodicity in the ROSAT data. Finding the same 
periodicity in the contemporary data of EGRET, as well 
as in old archival COS-B and SAS-2 data, yielded the 
value of the period derivative and thus of the object's en- 
ergetics. The discovery of the proper motion of the optical 
counterpart confirmed the optical identification and pro- 
vided the absolute position of Geminga to within 1" . Next 
came the measure of the source parallactic displacement, 
yielding a precise measure of its distance. Our knowledge 
of the Geminga pulsar is now good enough to warrant, 
" honoris causa" , inclusion in the radio pulsar catalogue 
Taylor, Manchester and Lyne (1993). In it, the radio silent 
Geminga stands out for the remarkable accuracy achieved 
in the measure of its parameters. For example, the period 
derivative of Geminga is known with an accuracy greater 
than that of the Crab. Such an accuracy is mainly due 
to the very stable behaviour of this 310^y old neutron 
star, which does not seem to be affected by glitches. In- 
deed, Mattox, Halpern and Caraveo (1996) claim that, 
during the first 3 years of coverage with EGRET, every 
pulsar revolution can be accounted for. This means that 
the 7-ray photons, collected in week-long observing peri- 
ods taken several month apart over a span of years, can be 
aligned in phase to form the well-defined, spiky light curve 
seen in high energy 7-rays. Had Geminga been a glitching 
pulsar, like Vela or Crab, it would have been impossible to 
obtain a satisfactory (and accurate) long-term solution. 
Moreover, the growing high energy coverage offers the pos- 
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sibility to refine the knowledge of this source's timing pa- 
rameters, inchiding the second period derivative and, thus, 
the neutron star bralfing index. However, to take full ad- 
vantage of the potential offered by 7-ray astronomy, very 
accurate values of the absolute source coordinates are re- 
quired to determine the barycentric arrival time of each 
photon. Any error in the source coordinates affects the 
accuracy of the corrected arrival times and thus hampers 
the global accuracy of the procedure. Mattox, Halpern 
and Caraveo (1996) have shown that the uncertainty in 
the source absolute positioning can induce an error in 
the barycentric correction up to 2.3 dg msec, where Sg 
is the source positional uncertainty (in arcsec) projected 
on the plane of the ecliptic. Thus, with a sensitivity to 
phase errors of < 10~^ over a period of 237 msec, the 
presently available positional accuracy of ^ 1" is the limit- 
ing factor for Geminga's photon timing analysis. As men- 
tioned before, high-energy 7-ray photons have been col- 
lected unevenly over the years: during 1972-1973 by the 
SAS-2 satellite, during 1975-1982 by the COS-B mission 
and again from 1991 by the Compton Gamma-Ray Ob- 
servatory. To lock in phase > 20 years of 7-ray data an 
improvement in positional accuracy of at least one order of 
magnitude is required. To achieve this challenging goal we 
have exploited the angular resolution of HST in conjunc- 
tion with the accuracy of the Hipparcos reference frame 
(ESA, 1997). 

2. From Hipparcos to HST 

However, to tie Geminga into the Hipparcos system, one 
has to overcome the 16-18 magnitude gap between the 
bright stars used as a primary reference and our niy = 25.5 
target, an impossible dynamical range for a single instru- 
ment. Moreover, one has to link astrometric images cov- 
ering a field of view of ~1 square degree to those obtained 
by the WFPC2 of HST, 14,000 times smaller in area. 
The principle of our method, described in detail by Car- 
aveo et al. (1997), is to cover the field of interest with 
images of increasingly smaller field of view and deeper 
limiting magnitude. The transfer of the reference system 
from one step to the next is based on sets of stars of inter- 
mediate magnitude. A similar procedure has been success- 
fully applied (e.g., Chiumiento et al, 1991; and Zacharias 
et al, 1995) for the determination of absolute positions of 
optical counterparts of extragalactic radio-sources. How- 
ever, these cases were less critical, as the differences in 
magnitude were only 8-10 and two or three intermediate 
steps were usually sufficient to bridge the gap between 
the reference stars and the targets. We note that the FK4 
coordinates of the 16 mag Crab pulsar were obtained by 
MacNamara (1971) with a two step procedure. 
In the case of the 9 magnitude fainter Geminga, we had 
to use a new five-step procedure starting with astromet- 
ric plates taken at the Osservatorio Astronomico Torino 
(OATo), and ending with images taken with the Planetary 



Camera 2 on HST (see Figure 1 for a comparison of the 

two images). The observational material actually used in 
this work is listed in Tablel. To calibrate the astrometric 
plates we have used the best optical reference frame in 
the sky which is now provided by the Hipparcos and Ty- 
cho Catalogues. These catalogues have been constructed 
in such a way that the Hipparcos reference frame coin- 
cides, to within limits set by observational uncertainties, 
with the International Celestial Reference System (ICRS). 
The latter system is defined by the adopted positions of 
several hundred extragalactic radio sources. It supersedes, 
although it is consistent with, the optical reference frame 
defined by the FK5 catalogue, which was formally based 
on the mean equator and dynamical equinox of J2000. 
While the Hipparcos Catalogue lists 118,000 stars as faint 
as m„ = 12.4 whose positions and annual proper motions 
are known to typically 1 mas at epoch J1991.25, the Tycho 
Catalogue contains more than one million objects, the as- 
trometric parameters of which have been measured with a 
median precision of 25 mas. With about 3 stars per square 
degree the Hipparcos Catalogue, although of vastly su- 
perior precision, is not suitable for our purpose, which 
can be better fulfilled by the relatively less accurate but 
denser Tycho Catalogue. However, since the proper mo- 
tions listed in the Tycho Catalogue are accurate only to 
approximately 25 mas/yr, we have propagated the Tycho 
positions (epoch J1991.25) to the OATo 38 cm refractor 
plate epoch (J1984.19) using the proper motions provided 
by the PPM catalogue (Roser and Bastian, 1993). This is 
an all-sky list of 378,910 stars referred to the FK5 system 
and provides proper motions with a typical accuracy of 
4 mas/yr. Consequently, we used a subset of the Tycho 
objects also listed in the PPM catalogue as our primary 
reference frame. On our first astrometric plate (left panel 
in Fig.l), 19 Tycho stars in common with the PPM define 
the primary reference, while a grid of 17 fainter stars is 
used as the secondary reference frame. This grid of fainter 
stars acts as the "primary reference frame" on the plates 
taken with the 105 cm astrometric reflector, defined as 
step 2 in Table 1. The procedure is then repeated to trans- 
fer the reference frame from step 2 to step 3, from 3 to 4, 
and from 4 to 5, ending with a grid of reference stars suit- 
able for the calibration of the WFPC2 frame (right panel 
of Fig.l). The optical and mechanical design of the OATo 
astrometric telescopes minimizes optical distorsions over 
fields of view even larger than the present ones. Indeed, 
no plate modeling residual effects (see i.e. Chiumiento et 
al, 1991; Lattanzi et al, 1991) can be seen in our data. 
The geometric distorsions in HST/PC data have been 
corrected following the standard procedure described in 
Holtzmann et al (1995). The accuracy of our procedure 
can be evaluated from the; values of er, the star centering 
errors which are typical of each telescope/detector combi- 
nation, and etr, the errors due to the least square proce- 
dure used to generate the secondary stars grid. These are 
also given Table 1. 
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Table 1. The stars in the secondary grid of step n (5th column) become the primary grid of step n + 1 (4th column). Since 
some of the secondary grid stars measured at step n could not be satisfactorily measured on the image of step n+1 (owing to, 
e.g., saturation, duplicity, etc.) the numbers in the two columns can be different. For each step we also give the centering (er) 
errors as well as the uncertainties arising from the frame transfer procedure (etr). The values for step 5 are educated guesses 
and not direct estimates, as for the other data. The limiting magnitude for each frame is listed in the last column. The value 
given for step refers to the stack of 10 frames of 15 minutes, each of them providing a limiting magnitude of 24-. 5 



step 
[N. images] 


telescope /detector 


field of view 
[epoch] 


primary grid 


secondary grid 




etr 


mag 


step 1 


OATo 38cm refractor 


70' X 70' 


19 Tycho/PPM 


17 


0".102 


0" .044 


14 


2 


plate, 30" /mm 


[f 984. 19] 












step 2 


OATo f 05cm reflector 


30' X 30' 


16 


28 


.061 


.020 


17 


2 


plate, 20".7/mm 


[1984.19] 












step 3 


OATo f 05cm reflector 


9' X 10' 


26 


21 


.026 


.011 


19.5 


2 


CCD, 0".48/px 


[1996.13] 












step 4 


ESO NTT 3.5m, SUSl 


2.5' X 2.5' 


16 


10 


.015 


.008 


26 


1 


CCD 0".13/px 


[1992.86] 












step 5 


HST 2.4m, PC2 


35" X 35" 


10 


Geminga 


[0 .005] 


[0 .003] 


26 


1 


CCD, 0".046/px 


[1995.21] 













20 arcmin 



Right Ascension 




Fig. 1. Comparison between the first and the last step of our procedure. Left: schematic representation of the 70' x 70' astrometric 
plate taken at the OATo. Hipparcos stars are marked with large filled ctrcles, Tycho stars with medium circles, while the small 
circles in the central region of the plate identify the stars used as secondary reference frame. The diamond in the center gives 
the position and actual dimension of the HST image, which represents the last step of our chain. Right: 35" x 35" Planetary 
Camera image of the field of Ceminga obtained with a 4400 sec exposure through filter 555, roughly equivalent to V. Ceminga 
is shown inside a white circle. 



The resulting position of Geminga at epoch=1995.21 

is 

a = 6''33'"54.1530", S = 17°46'12.909" ; these coordinates 
are in the Hipparcos reference frame and, therefore, in the 
ICRS. 

To this value we attach the error ^e|/sT + ^'sys + ^Tycho 
, where 

'^HST is the centering error of the HST image, which 



is conservatively taken as 1/10 pixel ( 0."005); 
'esys results from the combination of the errors induced by 
the least square adjustments used, at each step, to gener- 
ate the lists of secondary stellar grids Egys (Eichhorn and 
Williams, 1963). Following (Lattanzi et al, 1997), for m 
stars with average centering error Cr, such transformation 
errors can be written as etr ~ x Cr/^/rn, where 3 is 
the number of free parameters in the linear plate-to-field 
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transformation procedure. Table 1 gives the centering er- 
rors Cr and the frame transfer ones etr for each step of 
our procedure. Adding in quadrature the etrS one obtains 
Csys 0."050. Since this procedure has been independently 
performed for the two images available for steps 1 through 
3, a reduction of is to be applied to the single solution 

'^Tycho measures the precision with which wc can regis- 
ter the ensamble of 17 stars, defining the secondary ref- 
erence frame on the refractor image, relative to that re- 
alized by the 19 Tycho/PPM stars used as a primary 
reference on the refractor image. In our case tTycho ~ 
X (JTycho/V^, where arycho — 0."032 is the mean er- 
ror (per coordinate) at the epoch of the refractor plate ac- 
counting for both the average errors in the Tycho positions 
( 15 mas for our sample) and the uncertainty on PPM 
proper motions propagated over 7 years. Therefore, the fi- 
nal error to be attached to the ICRS position of Geminga 
is 40 mas per coordinate. 

Thus, the combined use of HST and Hipparcos yielded a 
25 fold improvement in the source absolute positioning. 

3. On the Use of an Accurate Position 

Among the pulsars seen in the optical or as high energy 
emitters, Gcminga's positional accuracy becomes the best 
so far, better than that of the 9 magnitudes brighter Crab 
pulsar (McNamara, 1971). This improved position, in con- 
junction with the HST measure of the proper motion and 
parallax (Caraveo et al., 1996), will allow the accurate 
calculation of the barycentric arrival time of all the avail- 
able 7-ray data and the locking in phase over more than 
20 yrs of data from three separate space missions. It will 
thus be possible to measure the period second derivative 
and, hence, the braking index n = vvjiJ^ of this neutron 
star. This quantity, expected to be 3 for magnetic dipole 
braking, has been only measured so far for young objects 
such as Crab (Lyne et al. 1988; n = 2.5 ± .01), PSR 0540- 
69 (Guiffes et al.l99; 2.8 ± .01), PSR 1509-58 (Kaspi et 
al.,1994; n = 2.0 ± .2), and for the slightly older PSR 
0833-45 (Lyne et al. 1996; n = 1.6 ± .3). 
While the results obtained for the three very young ob- 
jects are not too far from the expectations, the braking 
index recently measured for the ^ lO^yr Vela pulsar is 
definitely lower. Hence, the importance to exploit the sta- 
bility of the 3 10^ 2/ old Geminga to considerably enlarge 
the pulsar age sampled for braking index determination. 
The task is a challenging one : for a canonical braking 
index of 3, a i> of 2.710~^^s~^ is expected. This value is 
three order of magnitude smaller than that recently mea- 
sured for Vela and 6 order of magnitude smaller than that 
of the Crab. However, such tiny value of v is within reach 
of the SAS-2, COS-B and CGRO data. 
What makes Geminga suitable for the measurement of 
the frequency second derivative? If we order known pul- 
sars according to increasing values of v (expected for a 



braking index of 3), Geminga does not come out promi- 
nently. Not less than 40 pulsars have hypothetical v bigger 
than the value expected for Geminga. However, for all of 
them, these values are not measurable, because the errors 
on V and v are too big. What singles out Geminga is the 
possibility to reduce such errors by phasing together 20 
years worth of 7-ray data, now that the source positional 
accuracy is no longer a limiting factor. Thus, once again, 
Geminga appears to be in a special position amongst Iso- 
lated Neutron Stars. What matters most here is the source 
instrinsic stability, which made it worthwhile to devote a 
dedicated effort to the accurate measurement of its abso- 
lute position. 
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